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Abstract 

We propose a new observational test of top-down source models for the ultra-high- 
energy cosmic-rays (UHECRs), based on the simultaneous observation of two or 
more photons from the same Galactic hadron jet. We derive a general formula 
allowing one to calculate the probability of detecting such 'multiple events', for any 
particular top-down model, once the physical parameters of the associated hadron 
jets are known. We then apply our results to a generic top-down model involving 
the decay of a supermassive particle, and show that under reasonable assumptions 
the next-generation UHECR detectors would be able to detect multiple events on a 
timescale of a few years, depending on the mass of the top-down progenitor. Either 
the observation or the non-observation of such events will provide constraints on the 
UHECR top-down models and/or the physics of hadronization at ultra-high energy. 



1 Introduction 

Ultra-high-energy cosmic rays (UHECRs) are puzzling in respect of both their 
production and their propagation in the universe. On the one hand, even 
the most powerful astrophysical sites known to be able to accelerate parti- 
cles to very high energy seem to have difficulties to reach energies as high as 
3 x 10 20 eV (the highest reported UHECR energy so far[l]). On the other hand, 
even if they could, one would expect from the (presumably) extremely high 
rigidity of the observed UHECRs in the intergalactic medium that their di- 
rections of arrival in the Earth atmosphere roughly point towards the sources, 
which does not seem to be the case. Also, it had been expected that the 
UHECR flux above ~ 10 20 eV would be very much reduced due to the in- 
teraction of the UHE particles with the cosmological microwave background. 
This so-called GZK cutoff [2], however, does not seem to be present in the 
currently available data. 
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Although the UHECR sources are still essentially unknown, many models 
have been proposed, with various charms and problems [3]. They can be di- 
vided up into two classes: bottom-up models, in which particles initially at 
low (thermal) energy get accelerated by one or a series of astrophysical pro- 
cesses, and top-down models, in which each UHECR is directly produced, as 
a particle, at ultra-high energy, through the decay of a pre-existing super- 
massive particle or some exotic, high-energy physical process, e.g. involving 
the collapse or annihilation of topological defects. In this paper, we consider 
the so-called Galactic top-down models from a general point of view, merely 
assuming that the UHECR flux is dominated by sources in the Halo, with a 
density proportional to that of the dark matter. This solves the 'production 
problem' trivially (or more exactly shifts it to the problem of identifying the 
X-particles and explaining their production and decay rates) and provides a 
simple understanding of the absence of a GZK cutoff as well as the apparent 
isotropy of the UHECRs - at least until the statistics will be high enough for 
us to detect the dipole anisotropy due to the off-centered position of the solar 
system in the Galaxy (which should take at least three years of observation 
with the Pierre Auger Observatory (PAO) [5]). 

Many different models have been proposed, with X-particles of different types 
(either produced locally, notably through topological defect interactions, or in- 
herited from the big bang) and different masses (on the Planck scale, 10 28 eV, 
the GUT scale, 10 25 eV, or below) [3,6]. In this paper, we investigate a common 
consequence of a large class of Galactic top-down models, and propose an ob- 
servational test which could be accessible to the next generation of UHECR de- 
tectors, such as the PAO [7], the EUSO experiment [8] or the OWL/AirWatch 
project [9]. 



2 Multiple UHECR events: the basic idea 



Top-down scenarios involve the production of hadron jets in a way similar to 
what is observed in terrestrial accelerators, when an energetic quark- ant iquark 
pair (e.g. produced through e + e~ annihilation) hadronizes into a number of 
colourless hadrons through a QCD cascade. Let iV 7 be the number of photons 
in a jet (from neutral pion decay), and 0j et be the jet opening angle. Since 
gamma-rays propagate in straight lines away from the source, the average 
surface density of UHE photons at a distance D from the point where the 
X-particle decayed is a = N^/u^D 2 , where Uj et — n9f e t * s the jet solid angle. 
If a detector intersects such a jet, with a surface area S± orthogonal to the jet 
axis, it will see on average the following number of particles: 
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If the source is close enough, the number of photons in the jet high enough 
and the detector surface large enough, then fi will be larger than one and 
several UHECRs will be able to cross the detector at (almost exactly) the 
same time, from (almost exactly) the same direction. This is what we define 
to be a multiple event. It has an unambiguous experimental signature: two 
or more distinct showers developing simultaneously in the atmosphere, with 
almost perfectly parallel axes (within \fSZ/ D radians, which is much less than 
any conceivable experimental angular resolution). 1 

The number \i may be called the multiplicity of the X-particle decay event, 
or more exactly its potential multiplicity, as can be expected at Earth, since 
it is the average number of particles which can be observed simultaneously by 
the detector (assuming that it intersects the jet). For a given X-particle decay 
event, with a given /i, the actual multiplicity of the UHECR event as observed 
by the detector can only be predicted statistically. The probability, P(m, //), 
of observing an event with actual multiplicity m (integer) in a jet of potential 
multiplicity /i (real number) is given by the binomial law: 

(\ TO— 1 / \ A r 7 — TO 

where n/N 1 = Sj_/Sj Ct is the ratio of the detector's surface to the jet surface 
(see Eq. (1)), and thus the probability for a given particle in the jet to cross 
the detector. Note that V(m, ji) is actually the conditional probability of the 
multiple event, given the fact that one shower is observed, or if one prefers, 
the probability that a detected shower be accompanied bym-1 others. The 
probability of observing a multiple event with whatever multiplicity larger 
than two simply adds up to V(m > 2, /i) — 1 — (1 — fi/N^) N ' 1 ^ 1 wl - e~ M , for 
not too small values of N y . 

The basic idea behind top-down multiple events is thus that the UHECRs 
are not independent of one another, but appear in close groups released at 
the same time in a single X-particle decay event. If the groups are sufficiently 
tight, we should be able to detect several UHECRs at a time. In fact, top-down 
jets can be seen as genuine Galactic showers: in the same way as we detect 
single UHECR events by intercepting many secondary particles belonging to 
the same atmospheric shower, the use of very large detectors may allow us 
to detect single 'X-particle decay events' by intercepting several UHECRs 
belonging to the same Galactic shower. The detectability of multiple events 
thus comes down to the question: when we see one UHECR in a top-down jet, 
how close is the next one, compared to the detector's radius? 



Note that this is very different from the clustered events, sometimes referred to 
in the literature as multiplets, which correspond to independent UHECRs arriving 
from roughly the same direction in the sky, but at different times. 
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Fig. 1. Effective source density for various models of the dark-matter distribution in 
our Galaxy. Labels indicate the assumed value of R c for either isothermal or FRW 
models (see text). 

3 Timescale of multiple event detection 

The potential multiplicity, /x, of a UHECR event (Eq. 1) depends on two 
physical parameters related to the jet properties, iV 7 and 0j et , one astrophys- 
ical parameter, D, related to the source distribution, and one 'experimental' 
parameter, S±, related to the detector. Most X-particle decays will occur much 
too far from the solar system to give rise to multiple events. But if one assumes 
that the X-particles distribute over the Galactic halo in the same way as the 
dark-matter, one can estimate the probability that one of the many UHECR 
events that will be detected over a given period of observation corresponds to 
a small enough source distance. 

3. 1 The distribution of source distances 

The statistics of multiple events depend on that of source distances. For 
the dark-matter distribution in the Galaxy, we may consider either a simple 
isothermal halo model [11], where the density, udm depends on the galacto- 
centric distance, r, proportionally to l/(r 2 + R^), and the core radius R c is of 
the order of a few kiloparsecs, or an FRW model based on cold dark matter 
simulations, with udm l/[r(r + R c ) 2 ] [12]. In Fig. 1, we have plotted the cor- 
responding effective UHECR source density as a function of distance, for an 
observer located at the galactocentric radius of the Sun (taking into account 
the smaller effective contribution of more distant sources). 
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As can be seen, the effective source density is flat for low values of the source 
distance, which are those of interest to us because they give the highest prob- 
ability of observing multiple events. This result is nothing but the famous 
Olber's paradox, and it is independent of the actual dark matter distribution, 
provided it is not varying significantly on small scales. Moreover, the actual 
density profile of the dark matter halo appears not to affect significantly the 
normalization of the source density at small distances (except for unreasonable 
values of R c ). In the following, we adopt the value of 6 x 10~ 5 pc _1 (recall- 
ing that NRF models are currently preferred) and replace, for all practical 
purposes of the present study, the effective distribution of UHECR source 
distances by the following differential probability: 

dP(L>) = p(D)dD = -^-dD (0 < D < D ), (3) 
Do 

where Do = 1/(6 x 10~ 5 ) ~ 17 kpc is an effective radius beyond which no 
UHECR sources exist (that is, they contribute a negligible flux at Earth). 

Note that inhomogeneities in the dark matter distribution may in practice 
alter the probability of X-particle decay events at a given point of the Galaxy. 
A lower concentration of sources close to the Earth would decrease the chance 
of detecting multiple events, while a higher concentration would increase it. 
Lacking a precise knowledge of the small scale dark matter distribution, we 
can but assume that the Earth environment is not very different from the 
average. 



3.2 Multiple event probability 



For convenience, we shall rewrite the potential multiplicity of an individual 
UHECR event, given by Eq. (1), as: 

v( D ) = Hojjj, where /i = j^T^r ( 4 ) 

For any given model, the probability of detecting an event of actual multiplic- 
ity larger than n increases with the total number of UHECR events detected, 
N cvt , according to the simple law: 

V> n (N cvt ) = l-e-—, (5) 

where the constants N n are the characteristic numbers of events which have 
to be detected before it becomes reasonably probable (~ 63%) to detect an 
event of multiplicity n. This is a straightforward consequence of the statistical 
independence of X-particle decay events (see Appendix). 
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We show in the appendix that the characteristic event number for double event 
detection is given by : 

N2 = ^ 1/2 , (6) 
V 71 " 

and that the subsequent numbers deduce from N 2 by the following recursion 
relation (valid for detected multiplicities much smaller than the total jet mul- 
tiplicity, n<# 7 ): 

2n 

Nn+2 = 2—[ N n+li (7) 

so that, in particular, A3 = 2A 2 , N A = ±N 3 , N 5 = jjJV 3 , etc. 



3.3 Multiple event detection timescales 



In order to convert the above characteristic event numbers into multiple event 
detection timescales (for a given detector), we just need to calculate the 
UHECR detection rate. This depends on the total aperture, (in km 2 sr), 
and the duty cycle, 8 (in percent), of the detector. The UHECR detection rate 
above energy E is given by: 

A evt (> E) = $ CR (> E)xA d x5, (8) 

where <3>cr(> E) is the integral flux of UHECRs above energy E. From the 
AGASA and Fly's Eye experiments, a fair value of the UHECR flux at 10 20 eV 
is $cr(10 20 ) ~ 3x 10~ 40 cm -2 s -1 sr -1 eV -1 . The intergal flux obviously depends 
on the actual spectrum, which is virtually unknown above 10 20 eV. As for the 
spectrum below that energy, one should keep in mind that we are only inter- 
ested in the events which can be attributed to the top-down process under 
consideration, and which probably represent only a fraction of the total de- 
tected events between 10 19 and 10 20 eV. We shall note $xE th the corresponding 
integrated flux above the detector's threshold energy, E t h- 

We can now express the time evolution of the multiple event probabilities, by 
replacing A cvt by A cvt x t in Eq. (5): 

V> n (t) = 1 - e~t, where r n = N n /N evt . (9) 
Using Eqs. (6) and (8) and the expression for /z , Eq. (4), we find: 
9 iet (E th )D 2n 

T2= ^2^riS^' and T "^—i T " + >- (10) 

In practice, the average perpendicular surface of the detector, (S±), is related 
to the acceptance ^4d- If the detector's surface on the ground is Sd, and 6* max 
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is the maximum zenith angle visible by the detector, we have (S±) = ^S d x 
sin 2 6» max / (1 - cos6> max ), and A d = S d xtr sin 2 6> max . If max = 90°, (S±) = ±S d , 
and ^4d = irS d . Reporting into Eq. (10), we get: 

V2e iet (E th )D 



4 Numerical estimates for a toy jet model 

The remaining parameters necessary to calculate r 2 are the jet parameters, 
namely the photon multiplicity in the jet, iV 7 , or more precisely the number 
of photons as a function of energy, (dN^/dE)(E), and the jet opening angle, 
9 j Ct . Unfortunately, they are the most uncertain, because the detailed struc- 
ture of the hadron jets produced at the ultra-high energies of interest is not 
known, and one can only extrapolate from the semi-empirical models available 
at CERN energies, assuming that nothing dramatic occurs in physics at the 
intermediate energy scales. We shall not attempt here to describe QCD jet 
physics and theory, and refer the reader to the review by Bhattachargee and 
Sigl [3] of the various models, and to the book of Dokshitzer et al. [13], notably 
chapters 7 and 9, where the energy spectrum and multiplicity of the particles 
in a jet are discussed in detail, as well the collimation of both particles and 
energy. An interesting discussion of UHECR spectra in top-down models can 
also be found in [10]. 

4-1 The photon multiplicity in a jet 

Concerning the jet particle multiplicity, we show on Fig. 2a the typical spec- 
trum (multiplied by E 2 ) obtained with a modified leading- log approximation 
(MLLA) model, assuming that the X-particle mass is at the GUT scale, i.e. 
M x — 10 25 eV/c 2 (adapted from [3]). Also shown are the hadron jet spec- 
tra obtained with the Hill formula for X-particle masses of 10 25 eV/c and 
10 28 eV/c . In all cases, it is expected that most of the jet energy be dis- 
tributed among UHE particles. Following some previous works, we consider 
here MLLA spectra, which are found not to have a simple power-law behavior 
in the energy range of interest, namely around 10 20 eV, and to be steeper 
than often quoted for top-down scenarios - as would result from Hill's for- 
mula. It has been argued, however, that the MLLA spectra do not reproduce 
faithfully the fragmentation spectrum in the last energy decade or so, i.e. at 
energies close to the X-particle mass [10]. In this study, we shall only consider 
X-particle masses above 10 23 eV, so that the UHECRs of interest have energies 
well below m x . The MLLA spectral shape may thus represent a reasonable 
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Fig. 2. Left: Differential energy spectrum (multiplied by E 2 ) of the secondary pho- 
tons in a top-down jet, for various models (adapted from Bhattacharjee and Sigl 
[3]). Right: Integrated energy spectrum giving the total number of photons above 
energy E in the jet, for the same models. The slope of the approximate power-law 
(spectral index a) is indicated. 

approximation of the fragmentation spectrum around 10 20 eV. The normal- 
ization is probably more problematic, however, because it depends largely on 
the amount of energy carried out by the rarest, most energetic particles in the 
jet. Any other normalization than adopted below will lead to characteristic 
timescales for multiple event detection that scale according to Eq. (11), as 

In Fig. 2b, we have plotted the integrated hadron jet spectrum corresponding 
to the same cases as in Fig. 2a. If one tries to approximate the spectrum by 
a power-law dN/dE oc E~ x (i.e. N(> E) oc E~ a , with a — x — 1), the 
logarithmic slope is approximately constant for the Hill spectra, while it goes 
from 1.7 to 2.3 in the energy range of interest, with a value of ~ 2 at 10 20 eV, 
for the MLLA spectrum. In our 'toy jet model', we shall assume a mean 
hadron spectrum in E~ 2 between E ini = 10 19 eV and E sup = 10 22 eV: N(E) = 
T]Ej et /E 2 , where 77 is a numerical constant, so that N(> E) = rj(Ej Ct / E — 1) ~ 
TjEjet/E. A fit of the MLLA spectrum in Fig. 2 gives 77 ~ 0.10, so that about 
70% of the jet energy is in particles with energies above 10 19 eV. Only about 
one third of this energy, however, will be imparted to the photons, assuming 
that the total jet energy is divided up evenly into the three types of pions (tt°, 
7r + and 7r~), of which only the neutral ones decay into photons, and neglecting 
in a first approximation the contribution of nucleons. 

In conclusion, we shall adopt for the UHE photons the above E~ 2 spectrum 
with a value of r\ — 0.033. We shall also extrapolate the same spectrum to 
hadron jets generated by X-particles of lower mass, but with £j et , E- in i and 
E sup scaled linearly. This will allow us to explore UHECR progenitors with 
masses down to 10 23 eV. Assuming that the X particle decay events lead to the 
formation of two jets, so that E- ]et = M x c 2 /2, one finally obtains the following 
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approximate formula (only valid between and E sup ): 
4-2 The jet opening angle 

Coming now to the question of the jet opening angle, let us first note that a 
naive line of reasoning based on the Lorentz factor collimation effect cannot 
apply here. One might have been tempted to derive 6j et by claiming that an 
isotropic distribution of the jet particles in the rest frame of the parent quark 
would translate in the Galactic frame into a collimated distribution within a 
cone of opening angle 1/T 9 , where T q is the quark's Lorentz factor. However, 
such a collimation only applies to the decay products of real particles, not 
virtual ones. In the case considered here, the parent X-particle is at rest in 
the Galactic frame, and the jet particles are created out of the extremely 
intense field represented by a quark/anti-quark pair moving apart, not by the 
'decay' of one of its members. In a QCD jet, as it turns out, the hadronization 
process allows in principle large emission angles, i.e. large values of the particle 
momentum in a direction perpendicular to the jet axis, p±. 

A standard angular distribution is given by the following logarithmic law: 

™ I (13) 

dpl~ pl + A 2 Q cn [ } 

where the regularisation momentum, Aqqd, is a typical effective QCD scale, 
of the order of 300 MeV. Using this expression, and the fact that the emission 
angle of a particle is given by 9(p±) ~ P±/p\\ — p_i_c/E, one finds that about 
10% of the jet particles are found with p±_ < 30Aq C d, i-e. within 9j et ~ 9 x 
10 -11 (.E/10 20 eV)" 1 , with only a weak dependence on the X-particle mass 
(logarithmic in the large quantity Mxc 2 /Aqcd)- This represents a considerably 
weaker collimation than what would have been obtained from a Lorentz factor 
argument applied to a progenitor real quark, giving the opening angle 6- )et — 
l/r q ~ 3 x 10- 14 (M x c 2 /10 25 eV)- 1 . 

In practice, however, the energy is found to be better collimated than the 
multiplicity in QCD jets [13], and one expects the highest energy particles 
(which we are interested in) to be much better collimated than obtained from 
Eq. (13). In other words, the largest perpendicular momenta in the jet dis- 
tribution are attributed statistically more often to the lower-energy particles. 
Moreover, this behaviour is found to be amplified as the jet energy increases. 
This is important for our concern, because we are interested only in the par- 
ticles in the last few decades of the energy range (above a few 10 19 eV), and 
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not in all the (much more numerous) particles between this energy and the 
GeV range, which may fill a cone with a larger opening angle, but which shall 
not be detected as UHECRs anyway. 

Extrapolating the semi-empirical theory available at CERN energies towards 
the ultra-high energies of interest, one finds that for a quark jet, about 50% 
of the jet energy is found within #uhe ~ 2 x 10~ 12 radians of the jet axis 
(Dokshitzer, private communication). Considering this number as well as that 
obtained from Eq. (13), we shall arbitrarily adopt the following, hopefully 
conservative value for the UHECR jet opening angle: 

jet ~ 2 x 10" 11 . (14) 

The above estimates are admittedly arguable, and cannot be expected to hold 
for all models, but they may represent a reasonable description of the jets in 
the energy range of interest. Any other assumptions about d7V gamma /d.E and 
0j et (e.g. motivated by a detailed study of a particular top-down model) can 
be used in the following, in a straightforward replacement of ours. 

4-3 Observability of multiple events with the PAO and EUSO 

Let us now evaluate the characteristic timescale of double event observation, 
r 2 , by replacing the various model parameters by their numerical values in 
Eq. (11). The flux of top-down UHECRs above the detector's threshold energy, 
$> £ , is obtained consistently with the assumed hadronization spectrum: we 
normalize the E~ 2 spectrum to the quoted value of 3 x 10~ 40 cm _2 s _1 sr _1 eV _1 
at 10 20 eV, i.e. we assume that all the UHECRs at 10 20 eV have a top-down 
origin (and only a fraction of them below that energy). One thus obtains: 

$ cr(> E) ~ (10- 2 km- 2 yr- 1 sr- 1 ) x (j^—) , (15) 
from which it follows: 

f 2 1 1 x ( Sd \ ^ ( 6 \ 1 f-^U ^ ( J**-\ 
7-2-1 ' yT) V3000km 2 ; \100%J vl0 19 eW vl0 25 eW 

(16) 

In the case of the next generation UHECR observatories, the detection surface 
on the ground will be 3000km 2 for the PAO (one site), and 1.5 x 10 5 km 2 for 
EUSO. The detector's duty cycles are respectively 100% and 14%, and the 
energy thresholds are 10 19 eV for the PAO and 5 x 10 19 eV for EUSO. With 
these number, one finds, for an X-particle at the GUT scale (M x = 10 25 eV): 

r 2 (PAO) = 2.1 yr and r 2 (EUSO) = 0.48 yr, (17) 
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which are smaller than the observatories' lifetimes (15 and 3 years, respec- 
tively). The timescales for triple, quadruple and quintuple event detections 
are respectively 2 times, 2.67 and 3.2 times larger, as follows from Eq. (10). 

As indicated above, these timescales scale with the X-particle mass as M^ 2 , 
and with the actual number of UHE photons within the jets as N^ 2 . Even a 
drastic decrease in the photon multiplicity in the jet by an order of magni- 
tude would only increase the lifetimes by a factor \/T0 and keep the double 
event detection timescales smaller than the lifetime of each experiments. The 
sensitivity to 0j et is linear, however, and jet models with much larger opening 
angles than assumed here would make double event detection more problem- 
atic. We should also note that decay modes into more than two jets would 

1 /2 

lead to smaller photon multiplicities, so that r 2 would also scale as N- ' t . 



5 Strong upper limit on the double detection timescale 



Considering the above uncertainties, we shall now derive a model-independent 
limit for the detection of double events, based on the fact that every UHE 
photon in a top-down model comes from the decay of a neutral pion, and 
is therefore accompanied by a second photon within a very small angle, due 
to relativistic beaming. For a photon pair at 10 20 eV, say, the parent pion 
Lorentz factor is = 2E/m 7T c 2 ~ 1.4 x 10 12 , so that the opening angle of 
this minimal, two-particle jet is 9- sct ~ l/T n ~ 7 x 1CT 13 . Using the source 
distance distribution, Eq. (3), and averaging over the decay angle in the pion 
rest frame, one finds the probability distribution of the distance d between the 
two photons of a pair corresponding to a random UHECR event: 

p(d) = £lf for d«A)/rv (18) 

This allows us to estimate the probability that a detected UHE photon be 
accompanied by a second one within the range of the detector (of radius R d ): 

V 2 ~ JL T ^)± ~ 7 8 x ID- (^-) f-^U (19) 
2 "16 D "'- SXiU Uo 20 eVy v200W ' 1 ' 

From this, one can derive the minimum characteristic timescale for double 
event detection (independent of both 6- iCt and iV 7 ): 



Tr ~ N cvt V 2 ~ $ td (> E)^R\5V 2 ~ 23 yI V 14%J UokmJ ' (20) 

While this timescale may seem prohibitively long (17 years for the EUSO 
detector, of radius ~ 220 km), one should note the cubic dependence in R d : a 
detector only two times larger (i.e. on an orbit two times higher) would detect 
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double events from the minimal top-down jets on a timescale slightly above 2 
years. This model-independent upper limit may give confidence that multiple 
event detection should indeed be possible, if the UHECRs are produced in 
Galactic top-down jets containing not just two photons from an isolated tt° 
decay, but thousands of UHE photons. 



6 Conclusion 

In this paper, we have studied the possibility of observing multiple UHECR 
events from Galactic hadron jets, resulting from the decay of supermassive 
X-particles in the Halo. We have shown that, under reasonable assumptions 
about the jet properties, the next generation UHECR detectors should be able 
to detect a few double events, and possibly one or two triple and quadruple 
events, provided the UHE flux is dominated by top-down sources at about 
10 20 eV, and the mass of the X-particle progenitor is around the GUT scale. 

The main uncertainties in our calculations come from the jet model. Our as- 
sumptions about the photon spectrum, the jet multiplicity and the jet opening 
angle can only be considered as rough estimates, and other X-particle model 
may lead to different values. Nevertheless, we have derived a general frame- 
work for the study of multiple events probability, and given a way to calculate 
the relevant detection timescales for any X-particle model, once its physical 
parameters are specified. In particular, we have shown that the timescale r 2 is 
proportional to the jet opening angle, and inversely proportional to the square 
root of the photon multiplicity in the jets. 

Considering the jet uncertainties, we have also derived the double event de- 
tection timescale for the worst possible case: a two-particle jet consisting of 
the two photons produced by the decay of a neutral pion, as must be found in 
any Galactic hadron jet, whatever the model considered. This gives an upper 
limit on r 2 which would reduce to ~ 2 years for a detector twice as large as 
EUSO. Note that this is also independent of the X-particle mass, contrary to 
the timescales obtained by taking into account all the particles inside the jet. 
Note also that we have assumed homogeneously distributed photons inside 
the jets. In the case of a clumpy distribution, the probability of observing a 
multiple event can only be higher than what we have obtained here, because 
the photon density close to an arbitrary photon would then be higher (on 
average) than the mean photon density in the jet. 

Besides this new test of Galactic top-down models, three other observational 
signatures have already been proposed. First, top-down scenarios predict that 
photons should be the dominant component among UHECRs above 10 20 eV. 
As we have seen, this is very important for our study, because charged par- 
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tides of even ultra-high rigidity would be slightly deflected in the Galactic 
magnetic fields and lose the almost perfect collimation which multiple event 
require. Second, the dipole anisotropy due to the off-centered position of the 
Earth should eventually show up in the data, although this indirect evidence 
might not be fully discriminatory, since at least one bottom-up model has been 
proposed with the same characteristics [14]. Finally, the UHECR energy spec- 
trum should be characteristic of a hadronic fragmentation process, which may 
be quite different from the power laws usually expected from astrophysical 
acceleration processes. However, the exact shape of the spectrum is still hard 
to predict precisely in any of the top-down or bottom-up models, and it is not 
clear whether a measurement with reasonable error bars over two decades in 
energy at most (10 19 to 10 21 eV, say) can lead to definitive conclusions. 

By contrast, the detection of (were it only) one multiple event would provide a 
clear, direct evidence that a top-down scenario is involved, because it is highly 
improbable that two independent, bottom- up UHECRs arrive exactly at the 
same time from exactly the same direction. In this respect, we should recall 
that a previous study excluded (or more precisely found very improbable) the 
possibility that a heavy nucleus be photodisintegrated by the solar radiation 
field and give rise to a pair of showers from the lighter, daughter nuclei [15] 
(besides, such a pair would be perfectly correlated with the Sun's position). 

As shown above, the ideal detector to implement such a test and detect multi- 
ple events from Galactic-size hadron jets is a detector of very large acceptance, 
but not necessarily good angular and energy resolutions, which is usually one 
of the main challenges for UHECR detectors. One might therefore think about 
the interest of devising a detector made of a series of atmospheric fluorescence 
telescopes covering a huge surface on Earth, but with poor angular and energy 
resolutions in order to keep it economical. 

Finally, we note that while the jet parameters are a major cause of uncertainty 
in the calculation of multiple events detection timescales, this very fact may 
offer a possibility to constrain them (if the top-down origin of UHECRs were to 
be attested by this or another way). This may represent a unique opportunity 
to study the hadronization processes at energies many orders of magnitude 
above what can be reached in terrestrial accelerators. 
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A Appendix 



In this appendix, we derive the general expression of the probability of detect- 
ing a multiple event as a function of the total number of UHECRs detected, 
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which we shall note here N instead of N evt , for convenience. The probability 
of an event of multiplicity m > n will be denoted by V> n (N), and we shall 
prove that it can be written as in Eq. (5): 

V> n {N) = 1 - e-fc , (A.l) 
with the following values of the constants N n : 

N 2 = -Lio V2 , and N n+2 = -^-N n+1 , (A.2) 
Y 7T In — 1 

which hold for small values of n, compared to the jet multiplicity n <C iV 7 . 

Each useful X-particle decay event (i.e. giving rise to the detection of at least 
one UHECR), can be indexed by an integer k (1 < k < N), and is charac- 
terized by its distance, D k , to the detector. Its corresponding potential mul- 
tiplicity is given by Eq. (4): = nuD^/Df.. The statistics of multiple events 
detection will therefore be determined by the statistics of the source distances, 
Eq. (3): p(D) = 1/D , for < D < D . Combining both expressions, we find 
the probability distribution of the random variable \i: 

= ^oV 3/2 (a*o < A* < °o). (A.3) 



p(n) = p(D) 



dD 



dp, 



We shall first establish that, for one particular (useful) X-particle decay event 
with potential multiplicity p, the probability that it is not an event with actual 
multiplicity m > n writes, for n <^ N^: 

V(m>n,p) = (l + p + ^ + ...+ ^£^y) e"". (A.4) 

We proceed by recursion. For n — 2, the property reads V(m > 2,p) = 
1 — V(m > 2,/i) = e _M , and has been obtained in Sect. 2. Then: 

V(m > n + 1, p) = 1 - V(m >n + l,(i)= V(m >n,p) + V(n, p), (A.5) 

where V(n,p) was given in Eq. (2), and develops into: 

V{nfi) = (A6) 
F{n ^ } (n-l)\ A"™ -1 (l"^)"" 

and thus, for fi<iV r 

V(n,ri^j£-—e-^ (A.7) 
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which completes the proof. 

Now we consider the N events together, and note that the 'non-detection 
probabilities', V{m > n), simply multiply: 

N 

V(m>n,{fi k }) = Y[V(m>n,fi k ). (A.8) 
k=i 

Now the global multiple event probability, V> n (N), is the average value of 
V(m > n, {fi k }) = 1 - V(m > n, {//&}): 



V> n (N) = (l-l[V(m>n,fx k )j 

-i-(ni( 1+w . + ... + ^)e-»]) < A - 9) 

= 1 — (V(m > n, fi)^ , 
where the last equality follows from the statistical independence of the various 

This is indeed of the form announced in Eq. (A.l), provided that we define 
the characteristic numbers N n by: 

,~\ rv- (A.10) 

In ( P[m > n, [i) 

We thus have to calculate the average value of V(m > n, //): 

Vim > n, /i)) = (f>(m>n-l,n)) + ^-^yy (^"V") , (A.ll) 
where, with the probability law of Eq. (A. 3): 

fi n ~ 2 e^) = fi n ~ 2 e-^p{fi)dfi = ^- f°° //"- 2 - 3 / 2 e-M//. (A.12) 
' J no 2 J fj, 

Integrating by parts, one finds: 

V "V) = ^r 3 e""° + (n - 2 - I) (//""V") . (A. 13) 



We shall now use the fact that /i <l and limit the calculations to the lowest 

1 /2 

order (first order in /i ). We can thus drop the first term in the right-hand 
side of the above equation, as long as n > 4, and rewrite it as: 

(^-V) = (n - 2 - |) (^-V) , (n>4), (A.14) 
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from where it follows, using Eq. (A. 11) and writing I n = (V(m > n, fj,)j for 
simplicity, that 

n — 2 — - 

/„-/„_! = -^(4_!-4_ 2 ), (n>4). (A.15) 

n — 2 



This will allow us to calculate I n for all n > 4, once we know 7 2 and I3. 

-00 

IM) 



Starting with 7 2 = e ^p(fi)d/i and integrating by parts, one finds: 



J 2 = g-po _ ^2 ^ ^-i^-^ (A . 16) 

•//in 



and then, by changing the variable to u = fi 1 ^ 2 : 

h = e-« - 2^ /2 e - a du ~ e-"° - y^o 72 , (A.17) 

where we have used /i <^ 1 in the gaussian integral. Developing to first order 
in fi^ 2 , we thus obtain: 

J 2 ~l-v^/iJ /2 , (A. 18) 

and the first characteristic event number (for double event detection): 

N2 = f^-^ (A. 19) 

In I 2 y TT 

as announced in (A.2). 

Coming now to I 3 = J^(l + /u)e~ M p(/i)d/i, we have 

roc I , ("oo 

h = l2+ ^p(fi)dfi = I 2 + -fil /2 / ji-^e-Mji, (A.20) 

where we can replace the last term, using Eq. (A. 16): 

1/2 



(A.21) 



This allows us to write the characteristic number of events for the detection 
of triple events, A3 = — l/ln/3, as: 

A 3 = 4=/^o V2 = 2x4 (A.22) 



'7T 

From (A.18) and (A.21), we find 



W 2 = ^, (A.23) 
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which we can insert into (A. 15) to obtain: 



2 - 



and thus 



h = h + —f-ih g , 



3 ^Jix 3 



(A.24) 



(A.25) 



1/2 

More generally, we can show that, at the first order in /i , 

2(n-2) -1 
2(n-2) 



I n = 1 - a n y/7Tjjl /2 , with Q n = « n _ r 



(A.26) 



We have already found that this is true for n = 2 and n = 3. If we now assume 
that it is true for n — 2 and n — 1, we can use Eq. (A. 15) to calculate I n : 

I -I I ^-2) -3 



2(n-2) 

1/2 2n - 7 



1 - v^/V a n _i - 



= 1 - v / 7T/uJ /2 a n _ 1 

= 1 - y/TTfJ,l /2 a n -r 



4 K-1 - «n-2)v / 7T/ i /2 

2n - 7 2n - 7 2(n - 3) 



1 + 2n-4 2n-42(n-3)-l 
2fn-2) 



(A.27) 



2(n-2)-l' 
which is indeed the recursion relation (A.26). 



The characteristic number of events for a multiplicity larger than n is thus 
finally: 

JV„ = = — (A.28) 



In I n a n ypii ' 

and the recursion relation (A. 2) simply follows from that in (A.26). QED. 



18 



